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Resume : 

Un profit moyen turbulent est prescrit dans une conduite cylindrique, en adequation avec les observations experimentales. 
Nous considerons ensuite la nature des perturbations a cet ecoulement synthetique. Le calcul des croissances optimales 
predit deux types de structures, associees respectivement aux structures de proche-paroi et de grande echelle. Un ex- 
cellent accord quantitatif est trouve avec les resultats experimentaux quant a la longueur d'onde transversale. La reponse 
harmonique est egalement etudiee, et la croissance lineaire observee comparee a des simulations numeriques directes de 
turbulence forcee. Malgre de I'hypothese simple de type 'Eddy viscosity', cette approche lineaire predit ejficacement la 
croissance spectaculaire des modes de grande echelle au coeurde I' ecoulement. 

Abstract : 

A turbulent mean profile for pipe fiow is prescribed which closely matches experimental observations. The nature of 
perturbations superimposed upon this profile is then considered. Optimal growth calculations predict two distinct classes 
of structures, clearly associated with near-wall and large-scale structures. Quantitative correspondence of the spanwise 
wavelength of wall-structures with experimental observations is very good. The response to harmonic forcing is also 
considered, and the linear growth tested with direct numerical simulation of forced turbulence. Despite the very simple 
eddy viscosity assumption, this linear approach predicts well the surprisingly large growth of outer-scale modes in the 
bulk fiow. 
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1 Introduction 

The discovery of nonlinear solutions to the Navier-Stokes equations for Pipe and Couette geometries has led to 
significant progress of late. For pipe flow, these travelling wave solutions were discovered by adding a forcing to 
the Navier-Stokes equations to induce axially independent rolls. As the forcing is increased a three-dimensional 
wave instability arises. Three-dimensional states could then be traced back, using continuation methods, to the 
original equations upon reduction of the forcing HI 121 • Whilst successful at finding the first nonlinear states, 
the nature of this method leads to solutions that are only just self-sustained. As a consequence, the majority 
of known solutions (of the sinuous variety) were shown to be characteristic of the boundary between laminar 
and turbulent states. Initial states either side of this boundary develop into either turbulence or relaminarise 
ElSllllIB]. Only later were higher branches found (being also highly symmetric) with higher friction factors 
typical of turbulence |7 1. The highly organised structure of these states, however, does not appear to be typical 
of turbulence, nor do they display the strongly decelerated core of the mean flow. It is perhaps more intuitive, 
therefore, that one should begin with perturbations to this mean flow. 

2 Method 

Here we consider perturbations to a close approximation to the mean flow profile for pipe flow. The equations 
are normalised by the bulk flow speed Uh and the radius of the pipe R. Reynolds numbers are defined Re = 

2UbR/v and Rcr = UrRjv, where the wall velocity Ur = (^^ 5r.[/(r)) 2 |,.=ij is derived from the mean stress 
at the boundary. The mean profile U (r) is inverted from the 6,2;-averaged Navier-Stokes equations 

= -d,P + ;^ ( J + 5r] {yrd^fJ), B = -d,P, (1) 
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where the normaUsed effective viscosity is z^T(y) = 1 + E{y), and y = 1 — r is the distance from the pipe 
wall. The radially dependent eddy viscosity is prescribed ([8J, after [9]), 

Eiy) = ^ |l + ^^[2y - - + 2y'? 

We have updated the fitting parameters = 27 and k = 0.42 to be compatible with the observations of flOl 
and different scalings lead to the adjusted parameters R = Re / 2, B = 2 B. 

Perturbations, u, to the mean profile are first considered in the linear framework. The Navier-Stokes equa- 
tions with radially-dependent viscosity are linearised about U (y) ; note that the perturbations about U are also 
subject to the elevated effective viscosity. Rather than the usual progression to Orr-Sommerf eld-Squire form, 
eigenvalues and eigenvectors are found directly from the linearised primitive variable system with explicit so- 
lenoidal condition. Given the eigenfunctions and eigenvalues, using standard methods 1 1 1 1 one may calculate 
the optimal growth, 

G(a,m;t) = max ^['"'^^^^^ and Gmax = max (3) 

Wo \\Uq\\ t 

The wavelengths of the perturbations are = 211 /a (streamwise) and \q = 27r/m in azimuth (spanwise). 
Harmonic forcing is also considered in the following discussion. For a forcing of the form f{t) = /e"^*, a 
harmonic response u{t) = ue^'^* is expected. The optimal response is then 

R{a,m;Lo) = max — — and i?max = maxi?(L(j). (4) 
/ ll/ll 

The eigenfunction calculations were performed with up to = 200 points on r G [0, 1] using a Chebyshev 
collocation method. The number of (non-infinite) eigenvalues plus vectors obtained, being the number of de- 
grees of freedom, is 2 A^ — 3. Of these, 95% were kept for the optimal growth analysis. At the highest Re = 10®, 
for which Rer = 19200, the power spectral drop-off of the optimal mode was of 8 orders of magnitude (for 
both 'inner' and 'outer' modes, see below). 
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3 Results 

The structure of optimal modes in strongly sheared flows is well known to reflect the lift-up mechanism, 
whereby rolls aligned with the flow raise slow fluid from the boundary into the faster flow, leading to extended 
streaks. While shearing is certainly present in the turbulent flow, it is subject to a significantly larger effective 
viscosity, z/T(r). 

Optimal growths relative to the turbulent profile (figure [Til are significantly smaller than for the laminar case, 
although the greatest growth still occurs for the largest mode m = I. For the turbulent case, however, vt 
reduces to the laminar value as one approaches the wall, and growth of similar magnitude is also possible in 
the boundary region, provided that the length scale is not so small that diffusion again dominates. 
This secondary peak, for which rolls are located close to the wall, is shown in figure |2] to scale in 'inner' units, 
where in the absence of any other length scale, all units are inferred from the stress at the wall and the kinematic 




Fig. 1 - Left : Optimal growth at moderate Reynolds numbers. While the peak is for m = 1, it is clear that for 
larger a only modes of larger m persist. Right : At large Re a distinct peak appears at large m. 
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Fig. 2 - Left : An inner peak occurs for spanwise wavelength = 92. Right : The time s.t. G{T) = G, 




Fig. 3 - Structure of the optimal inner-mode (left) and its response (right). Each is normalised by the maximum 
of I It I . The initial condition consists of rolls, and the response is dominated by streaks. 



viscosity, r = r^u/ur, t = f^u/u^. It is clear from the figure that such scaling is appropriate at high Rcr, 
where the Gmax collapse with a peak for = 27r Rer/rn at 92. 

This peak in spanwise wavelength compares very favourably with experimental results. Using apparatus flush 
with the floor of the Utah desert [12] using smoke visualisation, a peak at was found at A+ = 93, for a Rcr 
as large as 1.5 x 10^. In [13J similar values (±10) were observed for Rcr in the range 1000 to 5000. Our 
predicted value A+ = 92 is in excellent agreement. In the present study, this secondary inner peak is observed 
to be distinct from the large-scale peak (m = 1) for Rcr of order 1000, and A^ being immediately close to the 
above value from the outset. Indeed, at lower Rcr, far before the peak becomes distinct for the a = case, it 
is clear that modes of finite streamwise extent select similar A^. 

In figure |3] the optimal initial condition and response are plotted for the case Re = 5300 (Rcr ~ 180) for 
m = 12 (A^ 98). The rolls are typical of those of the wall mode at much greater Re. 

The growth of an initial condition is all very well, but how easy is it to invoke such an initial condition ? The 
response to streamwise independent forcing is plotted over a range of harmonic frequencies, uj, in figurelH The 
steady forcing, lij = 0, is most effective at generating a response for this case. While this is not true for non-zero 
a, the peak response remained smaller than that for the axially independent case in all our calcuations. The 
code described in [14J has been modified to include radially dependent viscosity. The extreme difference in 
response to forcing different m modes was verified by timestepping relative to the prescribed mean profile. The 
structure of the optimal force /, however, is essentially identical to that of the optimal growth initial condition 
uq (see figure |3]l ; similarly for their respective responses. 

Unlike the transient growth, steady forcing provides a convenient platform for testing the model, whereby 
statistical measures averaged over much longer times may be accumulated. Consider the timestepped velocity 
field, u{t), obtained from direct numerical simulation of the original Navier-Stokes equations subject to a 

forcing /. Then Uf{t) = {u{t) ■ u)u is the component of u{t) in the direction of the the expected normalised 

response field u. The observed quantity may be compared with the expected response i?max- 

Figure \5\ shows the response of the turbulent flow Re = 5300, computed in a domain of length 10 radii 
(A^ f« 1800) at a resolution (60, ±64, ±64) before dealiasing. The calculations clearly demonstrate that large 
amplification may occur, as predicted by the model. Note also that for intermediate forcing, drag reduction is 
possible. 
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Fig. 4 - Response to harmonic forcing ; Rcr = 19200. For a = 0, steady forcing is largest {uj = 0). The outer 
mode (m = 1) is many orders more responsive for all moderate lo than the inner mode (m = 1314). 
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Response of a turbulent flow to forcing ; Rcr ~ 180, a = 0, m = 1 and ||/||^ = 10^^, 10^^, 10^^. 
IP = ^ J [ • I'^dV. The expected Rmax is 280. A similar order value is achieved for = 10"^ 

and decreases as the force increases. The amplitude of the induced flow is already as large as 5 to 10% of the 
mean flow for all cases. The forcing has a significant affect on the skin-friction coefficient C/ = 2{ur /Ub)^ , 
the horizontal line being the unforced average. 
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4 Conclusions 



The model has been shown to successfully predict large growth, to good quantitive agreement, despite the 
assumptions of linearity upon a turbulent state and the isotropic eddy viscosity assumption. It shows the huge 
response to forcing of large scale modes, and accurately predicts the peak spanwise wavelength, = 92, of 
the near- wall modes. 

Analogue studies for optimal growth in turbulent channels |[T5l[T6l . the boundary layer ifTTl and the Couette 
flow 1181 . all show close spanwise wavelengths for the wall mode. It is surprising, however, that the optimal 
growth calculations suggest small growth of the wall modes, despite clear observational evidence that such 
structures exist. 

The results from harmonic forcing, on the other hand, suggest that the large scale modes are easily excited. 
Such relative difficulty of forcing motion on small scales has consequences for control of turbulence. It clearly 
requires considerable effort to locally manipulate structures in the neighbourhood of the wall. It is possible 
that it is more straight forward to control indirectly via forcing of larger, or possibly even very large scales, 
taking advantage of the much greater linear response. Indeed, such a possiblity has been shown in principle for 
channel flow 1 19], and is verified here for a significantly larger ratio of 'large' (imposed roll) to 'small' (streak 
spacing) scales. The nature of suppression by a large-scale motion deserves further investigation. 
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